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Abstract 


Depending  on  the  initial  energies,  a proton  transfer  may 
proceed  either  via  a saddle- point  in  a potential  energy  sur- 
face or  by  crossing  from  the  reactants'  to  the  products' 
valley  before  the  saddle-point  is  reached.  In  the  second 
path  the  analogy  to  weak-  overlap  electron  transfers  is 
pointed  out.  The  present  study  is  intended  to  unify  pre- 
viously divergent  viewpoints,  by  showing  how  they  are 
special  cases  of  a more  general  picture.  Expressions  are 
obtained  for  the  reaction  rate  in  terms  of  the  properties  of 
the  potential  energy  surface  and  of  other  properties  of  the 
system,  using  a hydrogen  overvoltage  reaction  as  an 
example. 
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Introduction 


In  a proton  transfer,  for  example  in  the  hydrogen  overvoltage 
reaction  (1)  at  an  electrode  M, 

HsO+  + M(e)  - HaO  + H-M  (1) 

some  dynamical  effects  result  from  the  lightness  of  the  H-particle, 
and  several  questions  arise:  Can  a conventional  transition  state 
theory  be  used  to  calculate  the  reaction  rate,  with  a transition  state 
near  some  saddle- point  region  of  the  potential  energy  surface?  Does 
the  proton  transfer  occur  so  quickly  that,  as  in  weak- overlap  elec- 
tron transfers,  a "nonequilibrium"  solvent  dielectric  polarization 
arises? 
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Several  authors  assume  a conventional  transition  state  for 
reaction  (1)  and  also  an  equilibrium  solvent  dielectric  polarization. 
Others  assume,  instead,  concepts  analogous  to  those  used  in  a 
weak-  overlap  electron  transfer  reaction,  with  its  associated  non- 
equilibrium polarization.  Recently  a unified  treatment  was  out- 
lined,  and  a quantitative  description  is  given  in  the  present  paper. 

Potential  Energy  Surface  and  Reaction  Paths 

We  first  consider  the  potential  energy  surface  for  reaction  (1) 
as  a function  of  two  of  the  coordinates,  the  H-O-H  and  the  H-M  dis- 
tances,  using  the  usual  mass-weighted  skewed  axis  coordinates 
(Fig.  1).  There  are  also  the  bending  motions  of  the  O-H-M,  the 
stretching  motions  of  the  other  O-H  stretches,  and  the  coordinates  of 
the  solvent  environment.  The  actual  potential  energy  is  a function  of 
all  of  these  coordinates.  Path  a in  Fig.  1 is  a path  through  the 
saddle-point,  and  path  0 is  a path  at  any  fixed  O-M  distance. 

In  the  highly  exothermic  case,  caused  by  a very  favorable 
overpotential,  a schematic  diagram  of  the  surface  can  resemble  that 
in  Fig.  2.  Paths  a and  fi  are  again  drawn. 


Potential  energy  contour  plot  (schematic)  for  the 
three  center  reaction  H20-H-M  at  a fixed  value  of 
the  other  coordinates  and  at  a given  metal- 
solution  potential  difference,  for  the  case  where  the 
reaction  is  almost  thermoneutral  (symmetric).  The 
rotated  axes  are  scaled  H20-HM  and  H-M  distances 
(between  Q and  the  center  of  mass  of  HM  and  between 
H and  M).  9 The  configurations  along  the  dashed  line 
form  the  conventional  transition  state,  and  X denotes 
the  saddle-point  on  the  potential  energy  surface. 
Reaction  paths  a and  0 are  described  in  the  text. 


V n"M  "H,0-M 

Fig.  2.  Legend  as  in  Fig.  1,  but  the  reaction  is  now  highly 
exothermic.  The  dashed  line  denotes  the  conven- 
tional transition  state,  passing  through  the  saddle- 
point  of  the  potential  energy  surface. 


One  task  will  be  to  decide  whether  the  system  ever  reaches 
the  saddle-point  region  as  in  path  a in  Fig.  1 or  2,  or  whether,  either 
because  the  O-H  vibrational  energy  is  sufficiently  high  or  because 
tunneling  of  the  H from  the  O to  the  M may  be  large,  the  system 
moves  from  the  reactants'  valley  to  the  products’  valley  before  the 
saddle-point  region  is  reached,  as  in  path  0 in  Fig.  1 and  as  in  the 
0-path  at  a large  O-M  distance  in  Fig.  2. 

Paths  a and  0 can  be  competitive,  conceivably  one  (a)  pre- 
vailing in  very  exothermic  or  very  endothermic  conditions  and  the 
other  (0)  prevailing  under  more  nearly  thermoneutral  conditions.  A 
prescription  for  estimating  the  relative  importance  of  the  two  paths 
is  described  in  the  present  paper.  Path  0 is  assumed  in  Ref.  6 (but 
with  a surface  constructed  from  intersecting  parabolas)  and,  in 
effect,  path  a is  vised  in  the  usual  transition  state  theory. 

Energetics  for  Paths  a and  0 


We  shall  be  interested  in  introducing  a relatively  simple  for- 
malism which  allows  for  these  different  reaction  paths.  To  illustrate 
the  approach,  we  use,  for  simplicity,  a potential  surface  which 
treats  changes  in  the  O-H-M  potential  energy  due  to  changes  in  O-H 
and  H-M  distances  by  a bond  energy-bond  order  (BEBO)  method, 4,11 
and  which  treats  the  remaining  motions  largely  (though  not  neces- 
sarily entirely)  in  terms  of  their  dielectric  polarization  behavior.  In 
an  approximation  to  BEBO,  quite  adequate  where  tested,  the  BEBO 
surface  could  be  described  by  the  expression* ^ 

Ve  = ngAV0  4 AV*  [%! n nx  4 n^ln  nJAn  2 , (2) 

where  nt  and  n*  are  the  bond  orders  of  the  H20-H  and  H-M  bonds, 
respectively.  At  any  point  along  the  minimum  potential  energy  path 
it  is  assumed**’  *^  that  the  sum  n,  4 n2  remains  constant,  namely 
unity.  In  Fig.  3 we  have  Joined  by  a 0-path  any  pair  of  points  P and 
P*  lying  on  the  minimum  potential  energy  path,  but  in  the  reactants’ 
and  products’  valleys,  respectively.  Eq.  (2)  can  be  replaced  by  a 


more  elaborate  expression  without  changing  the  treatment  below. 


Fig.  3.  Legend  as  in  Figs.  1 and  2.  The  curve  denotes  the 

minimum  potential  energy  path,  which  proceeds  along 
the  reactants'  valley,  over  a saddle- point,  and  along 
the  products'  valley.  The  saddle-point  may  be  situ- 
ated as  in  Fig.  1 or  as  in  Fig.  2 or,  in  the  case  of  a 
very  endothermic  reaction,  in  the  products'  valley. 
Points  P and  P are  corresponding  pairs  points  wnich 
lie  at  the  intersection  of  a 0-path  with  the  minimum 
potential  energy  path. 


The  difference  of  bond  orders  of  the  H-M  bond  at  any  pair  of 
points  P and  P'  in  Fig.  3 will  be  denoted  by  An, 
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An  = n8(P/)  - aj(P)  . (3) 

Along  each  PP'  path  the  potential  energy  in  Figs.  1-3  is  given  by 
Eq.  (2)  (with  nt  + n,  = 1 at  the  points  P and  P')  and,  along  PP#,  by 
Eq.  (2)  using  bond  length-bond  order  relations. 


fj 


The  interaction  of  the  solvent  with  the  H2OH+-M  subsystem 
along  the  various  paths  is  also  to  be  included,  and,  for  brevity  and 
simplicity,  we  shall  do  so  classically.  (Various  high  frequency 
modes,  treated  structurally,  could  also  be  included  and  treated 
quantum  mechanically. ) We  consider  a particular  value  of  the 
orientation- vibration  polarization  of  the  remaining  coordinates  Pu(r) 
at  each  point  r of  the  solvent  medium. 
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An  expression  for  the  free  energy  of  solvation  of  the 
reactants  Gsol  for  any  given  Pu(r)  function  is13 

°soi  * -<»  - d;X/ ?■  sr  * ■ - s i-  f * * ; *•  p * «> 


where 


P ,,/D^ 

-vtf  op 


(5) 


Drt  and  D are  the  optical  and  static  dielectric  constants,  and  Dr  is 
is  the  electric  field  directly  due  to  the  charges.  Gg0l  varies  as  the 
point  P moves  along  the  minimum  potential  energy  curve  in  Fig.  3 
since  Dr  varies  with  position  along  that  line.  The  first  term  in  (4)  is 
the  solvation  term  when  there  is  no  orientation- vibration  polarization 
Pu,  the  second  term  is  a dipolar  interaction  of  Pu  with  the  charges 
in  a medium  of  dielectric  constant  Dop,  and  the  last  term  is  the 
orientation- vibration  polarization  energy  stored  up  in  the  polarized 
dielectric;  it  vanishes  when  D_  equals  D__ , as  does  P„. 
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When  the  system  in  Fig.  3 is  at  the  point  P the  solvation  free 
energy  is  that  of  the  products  G^ol  for  the  H20-H-M  configuration  P\ 
For  the  same  P„  it  is  given  by  the  same  expression  as  (4)  but  with  r 
subscripts  replaced  by  p's. 


In  the  vicinity  of  P in  Fig.  3 we  let  the  system  have  a local 
H20-H-M  protonic  vibrational  state  of  energy  E^,  a solvation  free 
energy  Gg0l,  and  an  electronic  energy  V*.  The  free  energy  of  the 
system  near  P,  Gr(P),  is  then  given  by 

Gr(P)  . Gg0j  * Ej  ♦ Vj  . (6) 

Ggoi  includes  the  electrode  potential  term  since  the  D in  (4)  includes 
fields  due  to  charges  on  the  electrode  and  in  the  ion  atmosphere 
(double  layer). 

Energetics  Alone  a fl-Path  (Proton  Jump) 


When  the  system  proceeds  along  a path  /3  in  Figs.  1-3,  from 


point  P to  point  P'  (Fig.  3),  it  does  so  by  a protonic  motion  so  rapid 
that  Pu  is  constant  along  the  PP'  path.  With  energy  conserved  during 
a protonic  jump  from  the  P to  P'  valley,  and  with  the  entropy  associ- 
ated with  the  orientation-vibration  polarization  Pu  also  unchanged 
during  the  transition  at  fixed  P„,  we  have 

Gr(P)  = GP(P')  . (7) 

Thus  far,  the  polarization  Pu  in  Eqs.  (4) -(7)  is  arbitrary.  As 
in  electron  transfer  reactions  we  choose  it  so  that  Gg0l  is  a minimum 
for  a system  at  the  point  P,  subject  to  the  constraint  imposed  by  (7) . 
Thereby,  one  finds 

6Gr  = 0 = - / (Dr  - 4 we  P)  • CP  dr  (8) 

6Gr  - = 0 = - f (Dr  - E>P)  • OP  dr  . (9) 

Multiplying  the  second  equation  by  a Lagrange  multiplier  m,  adding 
and  setting  the  coefficient  of  6PU  equal  to  zero,  as  the  most  general 
solution  of  (8)  and  (9),  one  finds  at  each  point  r in  the  medium  that 

4 ire  Pu/Dop  = (1  + m)  Dr  - m , (10) 

where  Dr  and  denote  the  electric  fields  for  systems  at  P and  at 
P/  in  Fig.  3,  directly  due  to  the  charges.  The  similarity  of  the  pro- 
cedure embodied  in  Eqs.  (4),  (8)-(10)  to  that  used^  for  the  transfer 
of  another  light  particle,  the  electron,  may  be  noted. 

Introduction  of  this  Pu  into  (4)  and  into  the  corresponding 
expression  for  G*?ol  yields  (11)  and  (14): 

C80l  * Glol  <«!>  + (11) 

where  the  first  term  is  the  equilibrium  solvation  for  a system  at 
point  P, 

Gsol(eq)  = "(1  * DsX)  / £r*  jf  d*/8* 


(12) 
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and  maA  is  the  fluctuation  term  due  to  Pu's  being  different  from 
its  equilibrium  value  in  Gg0l(eq).  A is  given  by  (13). 

A = f (DP  - Dr)  • (E>P  - Dr)  dr/8nc  . (13) 

The  difference  of  charge  distribution  on  the  right  hand  side  of  (13)  is 
expected  to  be  roughly  proportional  to  the  An  in  Eq.  (3),  and  so  to 
depend  on  the  point  P.  Similarly,  the  solvation  for  a system  at  point 
P'  in  Fig.  3 is 

<4*-  «5oi<*d  ♦«“♦»*»  • (H) 

The  value  of  m is  determined  from  (7),  (11),  and  (14): 

-(2m+  1)A  = GPQl(eq)  - G*ol(eq)  + E$  - Ej,  + V*?  - V*  (15) 

when  the  system  near  point  P in  a given  vibrational  state  v of 
energy  E*,  is  transformed  by  the  proton  jump  into  a system  near  P/ 
in  a proton  vibrational  state  v*  of  energy  E^r. 

Rate  Expression  for  the  g- Path 

We  denote  by  (e£)  the  probability  of  a reactive  v — V 
protonic  transition,  when  the  initial  translational  energy  along  the 
line  of  centers  in  Figs.  1-3  is  and  the  initial  protonic  vibrational 
energy  is  E°  The  transition  state  expression  for  the  reaction  rate 
is  given  by  * 5 


-(E°  + E?)/kT 


dEj  i/Q^  kT  (16) 


where 


f = kT  c~(AG8Ql*m*X)/kT  (2»iikT)  1 

^ h*  (2v#ikT)vk/h* 


= Z e 


-(ACjoi  + m*A)/kT 


(17) 
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AGgoi  is  016  increment  in  equilibrium  solvation  free  energy  on  going 
from  oo  to  P.  In  (17)  the  translational  partition  function  for  the  three 
translational  degrees  of  freedom  of  the  reactant,  and  that  for  a delo- 
calized reactant  on  the  surface  were  introduced.  (Corrections  for 
localized  adsorption  can  be  added  but  are  omitted  in  (16)-(17)  for 
brevity. ) Qy  is  the  protonic  vibrational  partition  function  for  the 
reactant  in  Figs.  1-3,  and  Z is  the  collision  frequency  (kT/2jr/i)  ^ 
for  collisions  with  unit  area  of  the  electrode. 


The  calculation  of  proceeds  as  follows:  One  first  locates 
the  point  P of  deepest  O-M  penetration,  namely  where  the  transla- 
tional energy  along  the  O-M  direction  vanishes,  L e. , where,  in  the 
reactants’  valley  we  have,  in  this  0-path  mechanism, 


Ey  and  AVg  refer  to  the  protonic  vibrational  energy  and  the  incre- 
ment in  potential  energy  (from  the  H30+  moving  from  »)  at  point  P 
in  Fig.  3.  In  calculating  E*  the  vibrational  quantum  number  v is 
taken  as  constant  (v)  within  the  reactants*  valley  (adiabatic  treatment 
for  the  H- vibration  in  the  valley).  Ey  differs  from  E°  only  because  of 
changes  in  cross-sectional  profile  (e.  g. , in  vibration  frequency) 
during  motion  along  that  valley  of  the  potential  energy  surface. 


If  the  value  of  Ey  at  P is  sufficiently  large  to  overcome  any 
Ve  barrier  from  reactants  to  products  along  the  0-path  at  that  P, 
the  corresponding  is  set  equal  to  unity.  Otherwise,  a tunneling 
calculation  for  is  used.  Recently,  a useful  calculation  which 
includes  tunneling  along  a 0-path  and,  where  necessary,  initially 
along  an  O-M  coordinate  was  given  in  Ref.  10  and  could  be  used  for 
the  present  purpose. 


The  barrier  along  this  0-path  is  modified  somewhat  by  the 
presence  of  the  Ggol  term  for  the  cited  Pu.  Thus,  changes  in  Ve  + 
along  this  0-path  serve  as  the  effective  barrier  to  proton  motion 


along  that  path  (Appendix). 


Thus  far,  the  electrode- solution  potential  difference  <p  has 
not  been  specifically  introduced.  It  is  implicitly  present  in  each 
Ggoi  terms.  By  evaluating  these  terms  one  obtains  (p.  As  in  elec- 
tron transfer  theory,*®  a "standard  potential"  <p°  for  reaction  (1)  in 
the  prevailing  medium  can  also  be  introduced  by  setting  the  free 
energy  of  activation  terms  AGr  and  ACjP  equal  at  that  <p.  One  can 
then  express  the  rate  of  the  forward  reaction  in  terms  of  the  dif- 
ference <p-<p°. 

Decision  as  to  Paths  g and  0 

If  for  any  point  P for  any  given  E°  and  E°  in  Eq.  (16)  the 
point  P is  closer  to  the  origin  than  the  saddle-point,  then  path  a 
will  dominate  rather  than  path  0 for  that  (E° , eJ)  pair.  If  point  P, 
in  the  case  of  Fig.  1,  is  quite  close  to  the  saddle- point,  the  An 
defined  by  (3)  becomes  very  small  and  so  the  term  Dr-D^  arising 
from  a difference  in  charge  distributions  at  P and  P'  also  becomes 
small,  and  so  does,  thereby,  the  X in  Eq.  (13).  Thus,  with  this 
approximate  vanishing  of  the  m2x  in  (17)  one  has  again  retrieved 
the  usual  transition  state  theory  result.  On  the  other  hand,  when  the 
barrier  along  a 0-path  can  more  easily  be  overcome,  either  by 
excess  vibrational  energy  in  E^  or  by  a sufficiently  large  value  of 
Kw# » one  stains  a 0-path  mechanism  rather  than  one  proceeding 
via  the  saddle-point  in  the  potential  energy  surface. 

The  above  remarks  also  serve  to  point  out  the  similarities 
and  differences  between  proton  transfers  and  weak- overlap  electron 
transfers.  The  weak- overlap  electron  transfer  proceeds  via  the 
counterpart  of  a 0-path,  and  has  the  m2X  terms  of  the  previous 
section.  The  proton  transfer  can  proceed  via  an  a -path,  wherein 
the  m2X  is  absent,  or  via  a 0-path,  depending  on  the  initial  con- 
ditions. The  remarks  made  earlier  on  the  hydrogen  overvoltage 
reaction  are  also  intended  to  apply  to  other  (e.  g. , homogeneous 


proton  transfers. 


Appendix:  Variation  in  D Along  a /3-P?*h 

Orf 

D varies  along  a /3-path,  since  the  electronic  structure  of  the 
system  varies  along  that  path.  D(r)  is  equal  to  (cf  Eq.  (5.  4)  of  Ref. 
13) 


D(r)  = -V,  J |*(rj)|,E 


ei 


£-£i 


n dr.  + C 
j ~J 


(Al) 


where  4>  is  the  electronic  wavefunction  for  any  nuclear  configuration, 
and  is  a function  of  the  coordinates  of  all  electrons  j of  the 
reactants,  and  the  sum  is  over  all  i,  L e. , over  all  electronic  and 
nuclear  charges  e.  in  positions  r ..  Thus,  this  D(r)  can  be  calcu- 
lated  from  an  electronic  structure  calculation.  The  C in  (Al)  is  the 
contribution  to  D arising  from  the  other  charges  on  the  electrode 
and  from  the  ion  atmosphere  (including  double  layer). 
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